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Abstract 

We consider cosmological implications of supersymmetric models with right- 
handed (s)neutrinos where the neutrino masses are purely Dirac-type. We pay 
particular attention to the case where gravitino is the lightest superparticle 
while one of the right-handed sneutrinos is next-to-the-lightest superparti- 
cle. We study constraints from big-bang nuleosynthesis and show that the 
constraints could be relaxed compared to the case without right-handed sneu- 
trinos. As a result, the gravitino-dark-matter scenario becomes viable with 
relatively large value of the gravitino mass. We also discuss constraints from 
the structure formation; in our model, the free-streaming length of the grav- 
itino dark matter may be as long as 0(1 Mpc), which is comparable to the 
present observational upper bound on the scale of free-streaming. 
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1 Introduction 



Existence of dark matter in our universe, which is strongly supported by a lot of re- 
cent cosmological observations [TJ [2j [3] , requires physics beyond the standard model. 
This is because there is no viable candidate for dark matter in the particle content 
of the standard model. Many possibilities of dark matter have been discussed in 
various frameworks of particle physics models so far [I]. Importantly, properties of 
the dark matter particle depend strongly on the particle physics model we consider. 

In this article, we adopt supersymmetry (SUSY) as new physics beyond the 
standard model. In such probably the most popular candidate for dark 

matter is thermally produced lightest neutralino which is usually assumed to be the 
lightest superparticle (LSP); in some part of the parameter space of the minimal 
supersymmetric standard model (MSSM), the relic density of the lightest neutralino 
well agrees with the present mass density of dark matter observed. However, as we 
discuss in the following, the lightest neutralino is not the only possibility of dark 
matter in supersymmetric models. 

If we try to build a supersymmetric model which accommodates with all theoret- 
ical and experimental requirements, we expect that there exist new exotic particles 
which are not superpartners of the standard model particles. For example, if local 
SUSY is realized in nature, gravitino tp^ which is the superpartner of the gravi- 
ton should exist. In addition, superpartners of right-handed neutrinos, which are 
strongly motivated to explain neutrino masses indicated by the neutrino oscillation 
experiments, may also exist. In particular, if the neutrino masses are Dirac-type, su- 
perpartners of the right-handed neutrinos are expected to be as light as other MSSM 
superparticles in the framework of gravity-mediated SUSY breaking. Importantly, 
one of those new particles may be the lightest superparticle and hence may be dark 
matter. In addition, existence of these exotic superparticles may significantly change 
the phenomenology of dark matter in supersymmetric models. 

In this paper, we consider the supersymmetric model in which the neutrino masses 
are Dirac-type and discuss cosmological implications of such a scenario. The possi- 
bility of the right-handed sneutrino LSP has already been discussed in Ref . [5] ; it has 
been pointed out that, if one of the right-handed sneutrinos is the LSP, the present 
relic density of the right-handed sneutrino may be as large as the dark matter 
density and hence the scenario of Dr dark matter can be realized. Here, we consider 
another case where the gravitino is the LSP and one of the right-handed sneutri- 
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nos is the next-to-the-lightest superparticle (NLSP). If the gravitino is the LSP, it 
may be a viable candidate for dark matter also in the case without the right-handed 
sneutrinos [6], [TJ, [8] . In such a case, however, stringent constraints on the scenario are 
obtained from the study of the gravitino production at the time of the reheating after 
inflation and also from the study of the big-bang nucleosynthesis (BBN) reactions. 
With the right-handed sneutrino NLSP, we reconsider cosmological constraints on 
the gravitino LSP scenario. We pay particular attention to the BBN constraints and 
also to the constraints from the structure formation of the universe. We will see that 
the BBN constraints could be significantly relaxed if there exists the right-handed 
sneutrino NLSP. In addition, we will also see that the free-streaming length of the 
gravitino dark matter may be a few Mpc, which is as long as the present sensitivity 
to the free-streaming length from observations. Thus, detailed understanding about 
the mechanism of structure formation has important impact on our scenario. 

The organization of this paper is as follows. In the next section, we introduce the 
model based on which our analysis will be performed. In Section [3], we discuss the 
BBN constraints on the gravitino LSP scenario with right-handed sneutrinos. We will 
see that the constraints could be significantly relaxed compared to the case without 
right-handed sneutrinos. Constraints from the structure formation will be discussed 
in Section HI Application of our discussion to the scenario where all the gravitino 
dark matter is produced by the decay of other superparticle will be discussed in 
Section Section is devoted to conclusions and discussion. 

2 Model Framework 

In this section, we summarize the model. As we mentioned, we consider the case 
where gravitino is the LSP while right-handed sneutrino is the NLSP. We assume 
that neutrino masses are purely Dirac-type, and the superpotential of the model is 
written as 

W = Wmssm + VuLH u O c R , (1) 

where JVmssm is the superpotential of the MSSM, L = (vlj&l) and H u = (H+,H° a ) 
are left-handed lepton doublet and up-type Higgs doublet, respectively. In this 
article, "hat" is used for superfields, while "tilde" is for superpartners. Generation 
indices are omitted for simplicity. In this model, neutrinos acquire their masses 
only through Yukawa interactions as m u = y u (H®) = y^vsm/3, where v ~ 174 
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GeV is the vacuum expectation value (VEV) of the standard model Higgs field and 
tan/? = (H®)/(H%). Thus, the neutrino Yukawa coupling is determined by the 
neutrino mass as 

2 \ 1/2 

\-13 ' 



y v sin (5 = 3.0 x 10~ ld x v — ? . (2) 

y H V2.8 x 10- 3 eV 2 / 

Mass squared differences among neutrinos have already been determined accurately 
by neutrino oscillation experiments [9j [10] , and are given by 

l Am l] atom ^ 2.8 x 10- 3 eV 2 , [Aral] solar * 7.9 x lO" 5 eV 2 . (3) 

In this article, we assume that the spectrum of neutrino masses is hierarchical, hence 
the largest neutrino Yukawa coupling is of the order of 10~ 13 . We use y v = 3.0 x 10~ 13 
for our numerical analysis. 

With right-handed (s)neutrinos, it is also necessary to introduce soft SUSY break- 
ing terms related to right-handed sneutrinos: right-handed sneutrino mass terms and 
tri-linear coupling terms called A^-terms. Soft SUSY breaking terms relevant to our 
analysis are 

£soft = -M\Vl - Ml R u* R i> R + (A V LHJ>% + h.c) , (4) 

where all breaking parameters, Mg, M„ R and A u , are defined at the electroweak 
(EW) scale. We parametrize A v by using the dimensionless constant a u as 

A v = a v y v M~ L . (5) 

We adopt gravity-mediated SUSY breaking scenario and, in such a case, a u is ex- 
pected to be 0(1). Though the A u -term induces the left-right mixing in the sneutrino 
mass matrix, the mixing is safely neglected in the calculation of mass eigenvalues 
due to the smallness of neutrino Yukawa coupling constants. Thus, the masses of 
sneutrinos are simply given by 

m 2 L = M 2 + icos(2/3)m|, ml R = M? R , (6) 

where mz — 91 GeV is the Z boson mass. In the following discussion, we assume 
that all the right-handed sneutrinos are degenerate in mass for simplicity. 

In this article, we consider the gravitino LSP scenario with right-handed sneu- 
trino NLSP. In such a case, the next-to-next-the-LSP (NNLSP) plays an important 
role in the thermal history of the universe. However, there are many possibilities 
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of the NNLSP, depending on the detail of SUSY breaking scenario. In our study, 
we concentrate on the case that the NNLSP is the lightest neutralino whose compo- 
sition is almost Bino. This situation is easily obtained if we consider the so-called 
constrained-MSSM type scenario It is not difficult to extend our discussion to 
the scenario with other NNLSP candidates. If Bino is the NNLSP, only the right- 
handed sneutrino which is the superpartner of the heaviest neutrino plays important 
roles in the Bino decay and other right-handed sneutrinos are hardly produced in 
this decay. Thus, mo R given in Eq. §6§ should be understood as the mass of the 
superpartner of the heaviest neutrino. 

3 Constraints from BBN 

It is well known that models with the gravitino LSP usually receive stringent con- 
straints from BBN. In these models, the lightest superparticle in the MSSM sector 
(which we call MSSM-LSP) has a long lifetime, sometimes much longer than one 
second due to Planck-scale suppressed interactions. Then, MSSM-LSP produced in 
the early universe decays into gravitino by emitting standard-model particles after 
BBN starts, which may spoil the success of BBN. In this section, we show that, when 
the right-handed sneutrino is the NLSP, the thermal history of the universe could 
be significantly altered, resulting in weaker constraints from BBN. 

Now, we consider the thermal history of the universe. In the early universe 
when the temperature is higher than the masses of MSSM particles, all the stan- 
dard model particles and their superpartners are in thermal equilibrium. Gravitino 
and right-handed sneutrinos are, however, never thermalized due to the weakness 
of their interactions. When the temperature becomes as low as (m^: mass 

of Bino-like neutralino), Bino-like neutralino B decouples from the thermal bath. 
Then, Bino-like neutralino decays into right-handed sneutrino as well as into grav- 
itino. After that, right-handed sneutrino decays into gravitino emitting right-handed 
neutrino. As can be easily understood, the decay of right-handed sneutrino is harm- 
less for the BBN scenario, because only right-handed neutrino is emitted. On the 
other hand, the decay of Bino-like neutralino into gravitino affects BBN. 

Constraints from BBN including hadronic decay modes are intensively studied 
in Refs.p^l IT3]: according to the studies, the BBN constraints give the upper bound 
on YxE v i S as a function of tx, where X stands for a long-lived but unstable particle, 
Yx = [nx I ' s\t^ Tx (with rix and s being the number density of X and the entropy 
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Figure 1: Two-body decay diagrams of Bino-like neutralino. 



density of the universe, respectively), E v i S is the mean energy of visible particles 
emitted in the X decay, and tx is the lifetime of X. We use the upper bound on 
YxEyis obtained in Ref.[12j. In our analysis, we adopt the line of 3^, (IT) in that work 
as the constraint from the p <-> n conversion. 



3.1 Decay of Bino-like neutralino 

First, we will take a closer look at the decay of the MSSM-LSP, which is Bino-like 
neutralino. Main decay modes are the following two-body decays shown in FigJTJ 
B — > vr&l, B — > tj}^, and B — ► ip^Z. Decay widths of these processes are given by 
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where 7723/2 is the gravitino mass and functions F and G are defined as 
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Here, gy is the U(l)y gauge coupling constant, M* ~ 2.4 x 10 is the reduced Planck 
mass, and 9w is the Weinberg angle. The lifetime of Bino-like neutralino is given by 
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Figure 2: Contour plot of the branching ratio Br(/3 — > ip^) (left figure) and the lifetime 
of B (right figure) on the (m 3 / 2 , rn§) plane: we take m^ R = 100 GeV, a v = 1, and 
mjj, = 1.5m^ in both figures. 

^ =T B- 2V B^u R D L + V B->^ + T B^Z' Where the factOT 2 in front ° f V B^u R v L 

comes from the contribution of the CP conjugate final state. 

The branching ratio of the mode B — > ip^ is shown in Fig[2] (left figure) on 
the (777.3/2, ma) plane, where we take m^ R = 100 GeV, a v = 1, and mo L = 1.5m§. 
Importantly, the decay mode B — > v^Pl competes with the mode B — > ^7 or it 
even dominates the total decay rate when the gravitino mass is larger than 0.1 GeV0. 
The lifetime of B is 10 2-3 seconds on most of the parameter region shown in Fig[5] 
(right figure). 

Without right-handed sneutrinos, B —>■ ip^/Z is the main decay mode, and sig- 
nificant amount of visible particles (including hadrons) are produced. As a result, 
the gravitino mass is strictly constrained as 777,3/2 ^0.1 GeV for < 1 second from 
BBN [13]. In our scenario with the right-handed sneutrino NLSP, however, less vis- 
ible particles are emitted, though the Bino-like neutralino is long-lived^. Therefore, 
constraints from BBN is expected to be relaxed. 

Next, we consider three- or four-body decay modes of B. Although branching 
ratios of these processes are much smaller than 1, they have impacts on the BBN 

x We have checked that z is about one order of magnitude smaller than rg_^, 7 on the 

parameter region of Figj2] 

2 Left-handed neutrinos injected by the decay might possibly change the abundance of 4 He [T3] . 
However, we have checked that the BBN constraints on the neutrino injection are much weaker 
than those on hadron injection from three- or four-body decay. 
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Figure 3: Three- or four-body decay diagrams of Bino-like neutralino. 



scenario. In particular, the hadronic decay modes give more severe constraints on the 
model than radiative ones via hadro-dissociation and p «-> n conversion processes. 

In our model, three- or four-body decay modes of the Bino-like neutralino are 
B — > ip^qq, B — > VRe^qq', and B — > VBpLqq shown in Figj3], where q and q denote 
quark and anti-quark, respectively. Since the hadronic branching ratio is compa- 
rable to radiative one in our scenario, the hadronic processes give the most severe 
constaints. Thus, we concentrate on hadro-dissociation and p <->• n conversion pro- 
cesses in order to derive BBN constraints. We calculate Shad^B-Evis as a function 
of Tg, where -Bhad is the hadronic branching ratio. In order to constrain the model 
quantitatively, we use the upper bound on Y x E vis obtained in Ref. [12ffi 

The product of the hadronic branching ratio and the visible energy, which is the 
mean energy of emitted hadrons from the three- and four-body decays, is given by 

B F ±- 

- D had- C/ vis -p 
1 B 



(9) 

where the factor 2 in first and second terms are from contributions of CP conjugate 

3 In some parameter regions where the Z boson in B — > ip^qq is off-shell, the photo-dissociation 
effect caused by B — ► ?/V7 might be comparable to the hadro-dissociation effect. In such a parameter 
region, however, B — > vrvl dominates in total decay as seen in FigO Therefore, the photo- 
dissociation effect is negligible, which is consistent with previous work [14] . 



^B^f, Z Bhad B vi 
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final states, -Bjf ad — 0.7 is the hadronic branching ratio of Z boson, (£v" is ) is the 

(z) 

averaged energy of hadrons emitted in each decay process, and E^ s is the energy of 
the Z boson, 

1/2 



E (Z) 
vis 




(10) 



3.2 Constraints 



In order to evaluate B^Y^E^, we have to determine the primordial abundance of 
Bino-like neutralino. The abundance depends highly on parameters in the MSSM 
sector such as masses of other superparticles. We use the following formula for the 
(would-be) density parameter of B [14] : 

2 



VLjjh 2 — C mo del X 0.1 



100 GeV 



ill) 



where the additional parameter C mo dei is introduced to take the model dependence 
into account: C mo dci ~ 1 for the neutralino in the bulk region, C mo dei ~ 0.1 for the 
co-annihilation or funnel region, and Cmodei ~ 10 for the pure Bino case without 
co-annihilation. Then, the yield of the Bino-like neutralino is given by 

YB = Cmodei x 3.6 x W^—^-^. (12) 

Our numerical results are shown in FigHJ where the BBN constraints are de- 
picted on the (m 3 / 2 , m B ) plane. We take Cmodei = 10, 1,0.1 in the left, middle, and 
right figures, respectively. Other parameters are the same as those used in Figj2j 
Shaded regions are ruled out by BBN. As shown in these figures, the constraints 
are drastically relaxed compared to those in models without right-handed sneutrinos 



As shown in the figures, new allowed region appears; for example, for C mo dei = 1? 
0.1 GeV < m 3 / 2 < 40 GeV. In that region, m B is bounded from above due to the 
BBN constraints from four-body decays, B — > VRe^qq' and B — > vrVlQQ- (Notice 
that the hadronic branching ratio -Bhad and mean energy E vis are enhanced when 
m B is large.) On the contrary, in the 0.01 GeV < ^0.1 GeV region, Bino- 

like neutralino decays mainly into the gravitino through the B —>■ ip^j process with 
the lifetime < 1 second. Since the decay occurs before BBN starts, it does not 
affect the BBN scenario. This situation also holds in the usual gravitino LSP scenario 
without right-handed sneutrinos, and the same allowed region can be seen in Ref . |T5] . 
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m 3/2 (GeV) m 3/2 (GeV) m 3/2 (GeV) 

Figure 4: Constraints from the BBN on the (m 3 / 2 , m §) plane: Parameters are chosen to 
be mc R = 100 GeV, a v = 1, and mp L = 1.5m^. We set C mo dei =10, 1, and 0.1 in the 
left, middle, and right figures, respectively. The middle-shaded regions are ruled out by 
BBN, while dark and light-shaded regions are excluded by the WMAP measurement and 
the structure formation of the universe, respectively. 

In the case of C mo d e i = 10(0.1), the constraints from BBN is more (less) stringent 
than the C mo dei = 1 case. As a result, the upper bound on m B becomes smaller 
(larger). Results also depend on left-right mixing angle 6p L -v R = \A u v/(m~ L — m| )| 
in the sneutrino mass matrix. As one can see in Eq.(7), branching ratio of the 
mode B — > ip^j/Z is more suppressed for larger 9p L .p R , which leads to less stringent 
constraints from BBN. On the contrary, when a u in Eq.(jSJ) is much smaller than 1, 
BBN constraints are severe and our constraints become close to those for the case 
without right-handed sneutrinos. We also find that the region m B ~ 100 GeV and 
0.1 GeV < m 3 / 2 < 1 GeV is excluded in the left and middle figures. In these resions, 
Bino is almost degenerate with right-handed sneutrino in mass. As a result, the 
process B — * urPl is kinematically suppressed and branching ratio of the process 
B — > ifipZ is enhanced. 

In addition to the BBN constraints, we also depict other cosmological bounds in 
FigJU the gravitino abundance originating in B must not exceed the value observed 
in the WMAP, Qnyih 2 ~ 0.105 [3J. Gravitino abundance from the decay is given by 
^3^2 = (^/wg) fig. Using Eq. ([TTj) . the cosmological constraint on the (m 3 / 2 ,wi B ) 
plane is obtained as m B m 3 / 2 < 10 4 GeV/C mo dei, which is shown as a dark-shaded 
region in FigJU This constraint gives the upper bound on m 3 / 2 - Another constraint, 
^3/2 < 0-4fiDM, is also depicted as a light-shaded region, which comes from the 
structure formation of the universe, which is discussed in the next section. 
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Figure 5: Diagram of the decay of vr. vr — > ^P^vr. 



4 Constraints from Structure Formation 

As shown in the previous section, larger value of 1713/2 is allowed compared to the 
case without right-handed sneutrinos. In the newly allowed parameter region, the 
MSSM-LSP decays mainly into right-handed sneutrino vr, and vr decays into the 
gravitino. Since the gravitino is produced with large velocity dissipation at the late 
universe, it behaves as a warm dark matter, and as a result, may affect the structure 
formation of the universe. In this section, we consider the constraints from the 
structure formation. 



4.1 Decay of right-handed sneutrino 



Once the gravitino is produced from the decay of vr, it is expected to freely stream in 
the universe and smooth out the (small scale) primordial density fluctuation. Since 
the lifetime of Bino is much shorter than that of right-handed scneutrino in our 
scenario, the free-streaming length of the gravitino is estimated as 



A 



FS 



dt 



"R 



a(t) 



2t E QU(t E Q) 




,(13) 



u 2 (t Dr ) 

where v(t) is the velocity, u(t) = p(t)/m%/2 with p(t) being the momentum of grav- 
itino, and a(t) is the cosmic scale factor. Here, the time of the matter-radiation 
equality and the lifetime of v R are denoted by £eq and t„ r , respectively. 

Right-handed sneutrino decays into the gravitino through the two-body decay 
process v R — > ip^UR shown in FigfSJ Notice that three- and four-body decay pro- 
cesses such as vr -> i)^ L Z, v R -> i) ll elW + ) v R -> ip^v L qq, and v R -> ip^qq' 
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are negligible, because these are strongly suppressed by the small neutrino mass 
compared to the two-body decay process. Thus, the decay width of i>r is given by 

r - t- 1 - 1 m ^ 

£*-^"* r ^ 4g7rM 2 m 2 /2 

The lifetime of right-handed sneutrino turns out to be 10 2 -10 8 seconds with mp R = 
100 GeV for 0.1 GeV < m 3/2 < 100 GeV. 

When m 3 / 2 is small enough compared to mo R , t eq ^> tj> r is satisfied and free- 
streaming length is approximately proportional to u(tp r )tI r . Since u{tz^)~ 1 and 

1/2 

t- r are both proportional to m 3 / 2 , Afs becomes independent of m 3 / 2 . With the use 
of the lifetime obtained in the above equation, the free-streaming length turns out 
to be Afs — 6 Mpc when ma R = 100 GeV unless m 3 / 2 is very close to my R . This 
fact indicates that the component of the dark matter (i.e., gravitino) from sneutrino 
decay acts as a warm dark matter (WDM). 

In addition to the sneutrino decay, gravitinos are also produced by the thermal 
scattering at the reheating epoch after inflation. The abundance of the gravitino 
from the scattering process is determined by the reheating temperature and the 
gravitino mass [Bl US] • Since the gravitino from the scattering is non-relativistic at 
the time of the structure formation, it acts as a cold dark matter (CDM). Thus, we 
have to consider the constraints from the structure formation of the universe on the 
WDM + CDM scenario. 



1 - 



m 



3/2 
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V R 



(14) 



4.2 Constraints 

Constraints from the structure formation on the WDM + CDM scenario are studied 
in recent works [TTJ [Tg] . According to these studies, it turns out that the matter 
power spectrum has a step-like decrease around the free-streaming scale of the WDM 
component, k ~ 27t/Afs- This is because only the power spectrum of the WDM com- 
ponent dumps at 27t/Afs- The power spectrum is estimated from the observations 
of the cosmic microwave background [H [3] , the red shift surveys of galaxies [2] , and 
so on. The WDM + CDM scenario is viable if the step-like decrease is within the 
uncertainty of the observed power spectrum. In this article, we adopt the following 
constraint: the power spectrum with the step-like decrease should be consistent with 
observational data p] at 95% confidence level. 

In our model, the energy density of dark matter is composed of two components, 
Pdm = P 3 / 2 + P 3 / 2 , where p%% and p^ 2 are the energy densities of gravitino produced 
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by the decay and by the thermal scattering processes, respectively. Introducing the 
fraction of WDM component /, we rewrite pdm as 

PDM = P3/2 + P3/2 = /PpureWDM + (1 ~ /)PpureCDMj (15) 

where p pur eWDM and p pur eCDM are the energy densities of pure WDM and CDM 
scenario, respectively, where n pure cDM^ 2 = ^pureWDM^ 2 — 0.1. We consider the 
adiabatic density fluctuation, then the power spectrum for the scale fc -1 is written 
as 

PMk) = [fPiH WDM (k) + (1 - /)P p 1 u / r 2 cC DM(^)] 2 • (16) 

In order to evaluate the magnitude of the step-like decrease, it is convenient to define 
the ratio of the CDM component in the total power spectrum at k > 27r/Aps : 

r ftM(^2,/A ) t _ /)2 (i7) 

-r pureCDMi^ ?C ^/Aps) 

The lower bound on r is obtained from the ratio of the lower and upper bounds on 
the observed power spectrum. In our scenario, Aps — 6 Mpc and hence k ~ 1 Mpc -1 . 
For such a wavelength, we obtain r > 0.35 at the 95% confidence level p], which 
leads to / < 0.4. In terms of the density parameter, it indicates ft^% < 0.4f2 DM , 
which gives the upper bound on the gravitino mass as m 3 / 2 < 40 GeV (4 GeV) for 
Cmodoi = 1 (10) (light-shaded regions in FigEJ). 

The constraint does not depend highly on the detail of the observational data. 
In fact, in other recent observations, it is claimed that the observational error on the 
power spectrum is about 15% [21 E], leading to the constraint as / < 0.2, which is 
of the same order of magnitude as the result above. 



5 Gravitino Dark Matter from Decay 

In our model, it is also possible to realize the scenario in which all the SuperWIMP 
(in our case, gravitino) dark matter is produced by the decay of other superparticle 
[7J [19]. In such a scenario, the SuperWIMP dark matter originates in MSSM-LSPc. 
It is well known that the relic abundance of the MSSM-LSP explains the observed 
dark matter abundance if the MSSM-LSP is stable. Therefore, even if the MSSM- 
LSP decays into the SuperWIMP at the late universe, the dark matter abundance 



4 Thus, it is implicitly assumed that the reheating temperature is low enough in order to suppress 
the gravitino production from the thermal scattering at the reheating epoch. 
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mg (TeV) mg (TeV) 

Figure 6: BBN constrains to the gravitino dark matter scenario on the (m^, ma R ) plane, 
where mp L = 1.5m^. The region m„ R > is irrelevant because we consider the case 
that ur is the NLSP. In the left figure, we set C mo d e i = 1, while C mo dei = 0.1 in the right 
figure. The constraints are shown as light-, middle-, and dark-shaded regions, for a v = 1, 
3, and 5, respectively. 



is still explained as far as the mass of the SuperWIMP is of the same order of that 
of the MSSM-LSP because of 

n _ ^SuperWIMP n ,-. ^ 

"SuperWIMP — "MSSM-LSP- {1°) 

m MSSM-LSP 

When the SuperWIMP is gravitino, the scenario receives stringent constraints 
from BBN [13], because not only gravitino but also visible particles are emitted in 
the decay of MSSM-LSP. In order to avoid the BBN constraints, one may assume that 
the lifetime of the MSSM-LSP is much shorter than one second or the MSSM-LSP is 
highly degenerate with the gravitino in mass. However, such possibilities often lead 
to fine-tunings of parameters in the model. On the other hand, in our model, the 
MSSM-LSP (i.e., Bino-like neutralino) decays mainly into the right-handed sneu- 
trino NLSP, then right-handed sneutrino decays into gravitino. The lifetime of B 
is shorter than that in model without right-handed sneutrinos, and the amount of 
visible particles emitted in the decay is suppressed. Furthermore, the decay of i> R 
into gravitino does not produce visible particles. Thus, the last decay is harmless for 
BBN, though the lifetime of sneutrino is much longer than one second. As a result, 
the constraints on the scenario are relaxed. 

In FigJHl we show the BBN constraints to the scenario on the Tnn R ) plane, 
where we take m^ L = 1.5m§. In the left (right) figure, we set C mo dei = 1 (0.1). Notice 
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that the gravitino mass is determined by Eqs.( |TTi) and ( jl~8i) with fisupcrWiMP — 0.1. 
The excluded regions from BBN are the light-, middle-, and dark-shaded regions 
for a u = 1, 3, and 5, respectively. In these figures, contours of the mass and the 
free-streaming length of the gravitino are also depicted. 

Since the gravitino is produced in the decay process, it may act as a warm 
dark matter. As a result, the scenario receives the constraint from the structure 
formation of the universe. From observation, the free-streaming length should be 
shorter than about 1 Mpc; with such a bound, we obtain the lower bound on right- 
handed sneutrino mass as mo R > 400 GeV. 

As shown in FigJHl when C mo d e i = 1, only the small region, ~ 600 GeV 
and mc, R ~ 400 GeV, is consistent with both of the constraints from BBN and the 
structure formation of the universe. The upper bound on comes from the BBN 
constraints on three- and four-body decays as discussed in previous sections. On the 
other hand, the upper bound on is due to those from the decay of B into the 
gravitino. When C mo doi = 0.1, the BBN constraints are relaxed, because the gravitino 
mass is larger than that in the C mo dei = 1 case and the decay rate into gravitino 
is suppressed. The constraints from the structure formation are also milder. As a 
result, a wide range of the parameter space is consistent with the constraints, which 
indicates that the scenario with gravitino dark matter from the MSSM-LSP can be 
naturally realized in our model if the annihilation cross section of the MSSM-LSP is 
large enough. 

Recently, it is pointed out that there are various discrepancies between numerical 
simulations for the structure formation based on the CDM scenario and the obser- 
vations of substructures in galaxies [201 l21j 122^ . Interestingly, a warm dark matter 



whose free-streaming length is slightly less than 1 Mpc may solve the discrepancies 
[T§1 12~4"l [23] . The free-streaming length as long as Aps ~ 1.0-0.4 Mpc is suggested 
as a solution to the missing satellites problem [20], which can be easily realized in 
our scenario. For cusp problem [21], however, parameter region which solves the 
problem receives noticeable change by analysis methods. For example, in |19j . it is 
claimed that the cusp and missing-satellites problem can be solved simultaneously 
when Afs is in appropriate range even if the lifetime of the decaying particle is as 
short as ~ 10 5 sec; then the cusp problem can be also solved in our scenario. On the 
other hand, in [241 . it is claimed that a simultaneous solution to both problems are 



5 There are also discussions that the observations using the gravitational lensing is quite consis- 
tent with not the WDM scenario but the CDM one [231- 
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hardly obtained unless the lifetime of the decaying particle is longer than ~ 10 10 sec; 
such a long lifetime cannot be realized in our scenario. The detailed analysis of the 
small-scale structure problems are out of the scope of this article and we leave them 
as future studies. 

6 Conclusions and Discussion 

In this paper, we have studied the cosmological implications of the gravitino LSP 
scenario with the right-handed sneutrino NLSP in the framework where neutrino 
masses are purely Dirac-type. In the case that MSSM-LSP is Bino-like neutralino, 
it mainly decays into the right-handed sneutrino with the lifetime ~ 10 2 -10 3 
seconds in the wide range of the parameter region. Though the MSSM-LSP is long- 
lived, no visible particles are produced in the leading process, thus constraints from 
BBN can be relaxed compared to the case without right-handed sneutrinos. With 
the quantitative analysis of the BBN constraints, we have found the new allowed 
region, 0.1 GeV < m 3 / 2 < 40 GeV, when mp R = 100 GeV. In this region, the BBN 
constraints give the upper bound on the Bino mass as < 200-400 GeV, which 
mainly comes from hadronic four-body decays. On the other hand, the upper bound 
on the gravitino mass is given by the constraints from the structure formation of 
the universe. In our scenario, some part of the gravitino is produced by the decay 
of right-handed sneutrino at the late universe. As a result, the gravitino freely 
streams in the universe and acts as a WDM. The gravitino is also produced from 
thermal scattering processes, which acts as a CDM. Taking the CDM contribution 
into account, we have considered the constraints on the WDM + CDM scenario 
from the observations of (small scale) structure formation, and finally found the 
upper bound on the gravitino mass. 

So far, we have concentrated on the case with Bino MSSM-LSP. Another well- 
motivated candidate for the MSSM-LSP is the lighter stau f x . With fj-NNLSP, 
the main decay mode of f\ is f\ — > vrV\%. We have also analyzed this case and 
derived upper bound on m^; we found < 400 GeV for m„ R = 100 GeV and 
1 GeV < m 3 /2 < 100 GeV. Here we take Uf L = 0.1, where U\l is left-handed stau 
component in lighter stau. In addition, the lifetime of lighter stau is 1-10 2 seconds 

6 If this process is kinematically forbidden, ti mainly decays as f\ — > ip^r and the right-handed 
sneutrino plays no significant role. Constrraints on such a scenario are already analyzed in [141126]. 
where the upper bound on the gravitino mass is given as 7713/2 < 10 GeV. 
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in the allowed region. Thus, the lighter stau may be seen as a long-lived charged 
track in future colliders. 

We have also discussed the possibility to realize the scenario where gravitino 
dark matter is produced from the decay of other superparticle. In the scenario, it 
is postulated that the gravitino dark matter originating in the Bino-like neutralino 
accounts for the total dark matter abundance. Considering the mass parameter 
region mo R = 400 GeV-1 TeV, we have found that the free-streaming length is Afs = 
1-0.3 Mpc, which allows to solve the small scale structure problems of galaxies. As in 
the case above, BBN constraints give the upper bound on the Bino mass as < 600 
GeV-1 TeV, which corresponds to m 3 / 2 < 80-250 GeV. As a result, all superparticle 
masses are within 100 GeV-1 TeV, thus the gravitino dark matter scenario in our 
framework seems to be natural. 
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